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O ; Abstract 

We propose, and constrain with the latest observational data, three model-independent 
parametrizations of the cosmic deceleration parameter q{z). They are well behaved and stay finite 
at all redshifts. We construct them by fixing the value of q at high redshift, q(z S> 1) = 1/2 (as 
demanded by cosmic structure formation), and at the far future, q(z = —1) = —1, and smoothly 

. interpolating q(z) between them. The fixed point at z = — 1 is not arbitrarily chosen; it read- 

>. 

IjTj , ily follows from the second law of thermodynamics. This fairly reduces the ample latitude in 

parameterizing q(z). 
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I. INTRODUCTION 



The deceleration parameter, denned as q = — 1 — (H/H 2 ), is a key quantity in describing 
the evolution of the homogeneous and isotropic universe. Its importance lies in the fact that 
it tells us the rate at which the Universe accelerates or decelerates its expansion. (Here and 
throughout H = a/a denotes the Hubble function, and a the scale factor of the Friedmann- 
Lemaitre- Robertson- Walker (FLRW) metric). Unfortunately, at present measurements of q 
suffer from non-small uncertainties that quickly grow with redshift (z = a -1 — 1 with ao = 1), 
though it is virtually certain that the Universe is accelerating nowadays, i.e., that qo < (the 
zero subscript means present time). Expressions of q(z) provided by cosmological models 
are of not much help either because none of them rests on sufficiently convincing theoretical 
grounds. It is to be hoped that things will eventually turn for the better when a successful 
theory of quantum gravity is in place, though it may well take a long while. 

Nevertheless, on the observational side the situation may improve comparatively soon 
given the variety and range of ongoing and planned major ambitious projects that involve 
bigger telescopes and advanced techniques -for a short review and a comprehensive list of 
references see §14 of [lj. In view of the above, it seems reasonable to propose parametrized 
expressions of q(z) based not in any specific cosmological model but on practical and em- 
pirical reasons that lessen their wide latitude. They can be useful while we patiently wait 
for a theoretically sound model backed by observation at all scales. 

Thus far, different parametrizations, such as q = qo + q x z, q = qo + qiz(l + z)^ 1 , 
q = qi + qiz{\ + z)- 2 , q = 1/2 + q x (l + z)~ 2 , q = 1/2 + {q x z + q 2 )(l + z)~ 2 , and 
more complex than these, have been considered in the literature to reconstruct q(z) from 
observational data (see e.g. [^-[sj]). However, the first parametrization is adequate for 
| z |<C 1 only and the others are unsuitable to predict the behavior of the deceleration 
parameter in the far future; in particular, they diverge as z — > — 1. Parametrizations whose 



intended range of validity includes the 



: ar future are necessarily more involved and usually 



contain three or four free parameters 



In this paper we propose three model independent parametrizations, with just two free 
parameters, valid from the matter dominated epoch (z ^> 1) onwards (i.e., up to z = — 1), 
based on practical and theoretical reasons and independent of any cosmological model. By 
construction they obey the asymptotic conditions, q(z ^> 1) = 1/2, q{z = —1) = —1, and 



a further condition, dq/dz > 0, which is valid at least when q — > — 1. The first condition 
expresses the conviction that at sufficiently high redshift the Universe was matter dominated 
(otherwise it would be very hard to account for the observed cosmic structures). At first 
sight, the other conditions are less compelling. As explained below, they are based on 
the second law of thermodynamics when account is made of the entropy of the apparent 
horizon. Usually one parametrizes a function in any specific interval by interpolating it 
between two given points (one at either end of the interval), modulo one first knows the 
value taken by the function at these two points. In actual fact, the parametrizations of 
q(z) proposed so far have just one fixed point: the asymptotic value at high redshift (g 
must converge to 1/2 when z ^> 1). The other, q , is not in reality a fixed point because 
the value of the deceleration parameter at z = is not very well known and therefore left 
free. The parametrizations proposed in this paper have two fixed points, one at the far past 
(z ^> 1), and other at the far future (z = —1). The second fixed point conforms to the 
thermodynamical constraints imposed by the second law. We believe this means a clear 
advantage over previous parametrizations of q{z), with just one fixed point. While it can 
be found in the literature parameterizations that also fix q at z = — 1 they do so arbitrarily, 
i.e., not grounded on sound physics. 

The aim of this paper is to propose three model-independent parametrizations of q{z), 
from z ^> lupto z = — 1, that comply with the second law of thermodynamics and constrain 
their two free parameters with recent observational data. As it turns out, all of them predict 
that the present stage of accelerated expansion will never slow down and are consistent with 
the spatially flat ACDM model. 

Section [Til considers the properties the deceleration parameter must observe assuming the 
Universe obeys the second law of thermodynamics. Put otherwise if, in the long run, it 
tends to an equilibrium state, i.e., of maximum entropy. Section ITUl presents the statistical 
tools, to be employed later, that make use of observational data from supernovae type la 
(557 data points), baryon acoustic oscillations combined with cosmic microwave background 
(BAO/CMB) (7 data points), and the history of the Hubble factor (24 data points). Section 
IIVI presents the three parametrizations, constrain them with the observational data, and 
discuss them. Finally, Section [V] summarizes our work and introduces some final remarks. 
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II. THERMODYNAMIC AL CONSTRAINTS ON q(z) 



As is well known, physical systems tend spontaneously to some equilibrium state com- 
patible with the constraints imposed on them. This summarizes the empirical basis of the 
second law of thermodynamics. Very briefly, this law establishes that isolated, macroscopic 
systems, evolve to the maximum entropy state consistent with their constraints 11] . As a 
consequence their entropy, S, cannot decrease at any time, i.e., dS > 0. Further, in the last 
phase of the evolution S has to be a convex function of the said variable, d 2 S(z — > — 1) < 0. 

Arguably, the entropy of Friedmann-Lemaitre-Robertson- Walker (FLRW) universes is 
dominated by the entropy of the causal horizon, at least at late times -see e.g. 12] . As 
causal horizon we shall take the apparent horizon, the marginally trapped surface with van- 
ishing expansion of radius f ^ = 1 / V H 2 + ka~ 2 [l3| , where k denotes the spatial curvature 
index. Interestingly enough, it has been shown that this horizon represents the appropri- 
ate thermodynamic boundary surface jl^ . Leaving aside possible quantum corrections its 
entropy results proportional to area of the latter (47rf^) [131 ]. 



S A oc A = 4vr (H 2 + k a 



(1) 



Therefore, so long as we can ignore the entropy within the horizon, the second law of 
thermodynamics imposes A! > 0, at any time, as well as A" < at late times -the prime 
meaning derivative with respect to the scale factor. Both conditions are to be fulfilled if the 



15| 



FLRW universe is to tend to thermodynamic equilibrium at late times 

Bearing in mind the definition of the deceleration parameter, for spatially flat (k = 0) 
FLRW universes we can write 



A' = 2A 



l + q 



and 



A" = 2A 



<1_ 

a 



l + q 



(2) 



The first equation implies q > — 1. Inspection of the second one reveals that when a — > oo 
the middle term in the square parenthesis dominates. Thereby, dq/da < in that limit. 
Thus, for the Universe to tend to thermodynamic equilibrium at late times we must have 
q — )■ — 1 and dq/dz > as z — > —1. 



III. STATISTICAL TOOLS 



This Section outlines the use of the observational data to fit the parametrizations of the 
deceleration parameter, q{z). Since the likelihood function is defined by C oc exp(— x 2 /2) 
the best fit to the data follows from minimizing the sum Xtotai = Xsn + Xbao/cmb + Xh- 
As detailed below, the best fit values of the parameters can be obtained by contrasting the 
proposed parametrizations with the empirical data mentioned above and minimizing the 
Xtotai by means of the Markov Chain Monte Carlo method. 



A. SN la 



We compare the theoretical distance modulus 



^th(zi) = 51og 10 + , (3) 

where [Aq = 42.38 — 51og 1Q /t, with the observed distance modulus ii bs{zi) of the 557 su- 



pernovae type la assembled in the Union2 compilation [16|. The latter data set is substan- 
tially richer than previous SN la compilations and presents other advantages; mainly, the 
refitting of all light curves with the SALT2 fitter and an upgraded control of systematic 
errors. In (j3]) dj, = (1 + z) J* E ^. p ^ is the Hubble-free luminosity distance, p = (q±, #2), and 
E(z; p) = H(z\ p)/H . Here and throughout, qi and q<i denote the free parameters occurring 
in the parametrizations presented below in Section IV; H is the Hubble constant and h its 
value in units of 100 km/s/Mpc. 

The x 2 from the 557 SN la is given by 

XsnW-^ ^ , W 

where the subscripts u th" and "06s" indicate the theoretical value (i.e., the value from the 
parametrization) and the observed value, respectively. As usual, the <j{ quantities stand for 
the la uncertainty associated to the ith data point. To eliminate the effect of the nuisance 
parameter /x n , which is independent of the data points and the data set, we follow the 
procedure of [f 7| . 



B. BAO and CMB 



Baryon acoustic oscillations can be traced to pressure waves at the recombination epoch 
generated by cosmological perturbations in the primeval baryon-photon plasma. They have 
been revealed by a distinct peak in the large scale correlation function measured from the 
luminous red galaxies sample of the Sloan Digital Sky Survey (SDSS): at z = 0.35 [181 ]. as 



well as in the two degree Field Galaxy Redshift Survey at z = 0.2 



peaks have been observed: at z = 0.278 (with the SDSS [20|), at z = 0.106 (in the six degree 



191 ] . More recently other 



2l|), and at z = 0.44, z = 0.60, and z = 0.73 (by the WiggleZ 



Field Galaxy Redshift Survey 
team j^J). 

From each peak the ratio of the comoving sound horizon r s (z) = c s (z)/ H(z) dz at de- 
coupling (z = z± ~ 1090) and at the drag epoch (z = Zd), the epoch at which the acoustic 
oscillations are frozen in, can be measured. Here c s is the speed of sound. Likewise, at each 
peak, a characteristic distance scale, the dilation scale 

d\(zBAO 



Dv(zbao) 



ZBAO 



H(z BAO ) 



(5) 



2 dz' 



is the comoving angular distance, can also be determined. To 



where d A (z) = J H{zl) 
compute the drag epoch redshift we use the formula (4) of Eisenstein and Hu in 23( and 
get z d ss 1020. 



Multiplying the ratio 



, d^zb\ ) ' ^ a ^ en f rom the BAO peaks by the acoustic scale 



we get d ^( z *) r 3 (z d ) ^ eacn rgcig^ift Q f the seven BAO data. Here we use the value for 1* 

to D v (zBAO) rs(z*) A 



derived from Wi 
302.09 ±0.76 



kinson microwave anisotropy probe (WMAP)7-years data, namely, I a 



241 ] . If we also use the value of the ratio of sound horizon at the drag epoch 



and at recombination (redshift z+), computed from the values reported in 



24 



r s (zd) 



1.045 ±0.015 we obtain the new estimator 



dA{z+) 



., shown in table [Tl as done in ^5|. Using 
this estimator, the dependence in the sound horizons at decoupling and the drag epoch is 
suppressed. Thus we just use the ratio between them, which is almost model independent. 
This follows because both redshifts are rather close and the sound horizon at decoupling 
and drag essentially depend on the fractional difference between the number of photons and 
baryons [25j. To obtain the x 2 f° r the combined BAO/CMB data we compute 



ZBAO 


0.106 


0.2 0.278 0.35 0.44 


0.6 


0.73 




0.336 ±0.015 


0.1905 ±0.0061 0.1394 ±0.0049 0.1097 ± 0.0036 0.0916 ± 0.0071 


0.0726 ± 0.034 


0.0592 ±0.0032 


Dv(znAn) 


£*a(z*) 


30.92 ± 1.45 


17.53 ±0.62 12.83 ±0.49 10.09 ± 0.36 8.43 ± 0.66 


6.68 ±0.33 


5.45 ±0.30 


°1 



Table I. Values of 



r a (z d ) 
Dv(zbao) 



(reported in [184221]) and the derived ratio jj 



dA(z*) 
v(zbao) ' 



where 



Xbao/cmb — X T C X X 



/ d A {z+) _ on qr>\ 
' LV(0.106) t)U - yZ » 



(7) 



X 



V 



dA(z*) _ 
D v (0.2) 
dA(z*) 
Dy(0.278) 

dA{z*) 
D v (0.35) 

dA(z*) 
D v (0.M) 

dA{z*) 
ZM0.6) 

dA{z*) 
D v (0.73) 



17.53 

- 12.83 

- 10.09 
-8.43 
-6.68 
-5.45 J 



and X T the transpose matrix. The elements of covariance matrix C are given by 



' f) d A (z±) ■ 
U D v (z) 



' f) d A(z*) ' 

U D v (z) 



dp k I \ dp k 
where the sum is over the estimators used fin our case 



r a {z d ) 
' D v {z b ao) 



, l A /n and 



(8) 



The 



elements of the original covariance matrices are Ci A — ffj, , Cr, 



(*d) 



e 2 



and C B ao s 



rs(z d ) ^BAUlJ — 

rs(z*) r a (zi<) 

r_A*A, Tll e 



9b ao i 9b ao j r ij where Qbao i stand for the errors associated with the estimator 
only non-zero off-diagonal correlation coefficients r^- are r^_ 02z= o.35 = 0.337, r z=0A4:Z=0l , 



19 



221 ] . Thus, the inverse covariance 



0.369 and r z=o ez=Q 73 = 0.438, and their symmetric 
matrix comes to be 

/ 0.492 -0.084 -0.126 -0.136 -0.025 -0.081 -0.088\ 

-0.084 3.362 -0.327 -2.397 -0.065 -0.209 -0.228 

-0.126 -0.327 4.429 -0.528 -0.098 -0.314 -0.342 

-0.136 -2.397 -0.528 9.712 -0.106 -0.338 -0.368 

-0.025 -0.065 -0.098 -0.106 2.798 -2.749 1.182 

-0.081 -0.209 -0.314 -0.338 -2.749 15.002 -7.294 

\-0.088 -0.228 -0.342 -0.368 1.182 -7.294 14.587 J 



C 1 
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C. History of the Hubble parameter 



The history of the Hubble parameter, H(z), is poorly constrained though, recently, some 
high precision measurements by Riess et al. at z — 0, obtained from the observation of 
240 Cepheid variables of rather similar periods and metallicities j^fjj], and Gaztanaga et al. 
at z — 0.24, 0.34, and 0.43 2t|], who used the BAO peak position as a standard ruler in 
the radial direction, have improved matters somewhat. We have employed these four data 
alongside 11 less precise data, in the redshift interval 0.1 < z < 1.8, from Simon et al. 



28l | and Stern et al. 29|, derived from the differential ages of passive-evolving galaxies and 



archival data. In addition we have included in our analysis 9 more recent correlated data 



from the WiggleZ survey 
The corresponding \ 2 



30]. 



X 2 H (P) 



15 

£ 

i=l 



[H th ( Zi ) - H obs ( Z j)f 



(9) 



where 



X 



H 



f H th (0M) 


- 69.4^ 


H th (0.l5) 


- 76.6 


H th (0.2h) 


- 75.3 


H th (0.35) 


- 78.3 


Hth(0A5) 


-87.3 


H th {0.55) 


-88.9 


H th (0.65) - 


- 101.4 


H th (0.75) 


-96.9 


\H th (0.85) - 


- 127.3y 



and C^ 1 is the inverse covariance matrix given in table 6 of 30]. 



IV. PARAMETRIZATIONS 



Here we propose and constrain three parametrizations of the deceleration parameter, 
valid from the matter dominated era up to z = —1. These fulfill: (i) q(z ^> 1) = 1/2 
(as demanded by cosmic structure formation), (ii) q{z = —1) = —1 and dq(z)/dz > 
when q(z) — > —1 as required by the thermodynamic arguments of above (the second law of 

8 



thermodynamics). In interpolating between [z = — 1, q = — 1) and (z > 1, q — 1/2) we 
introduce two free parameters q± and q2 and fit them to the observational sets of data by the 
method of last Section. Note that due to the scarcity of q(z) data and their big error bars, 
we do not constrain the parametrizations directly. We constrain instead the expressions for 
H(z) that arise from integrating them; namely, 

H(z) = H exp jjf [l + g (x)]dln(l + x)| , (10) 

which holds for all parametrizations. This has the advantage of a much bigger and robust 
statistics. Notice that the Hubble constant also enters this expression as a free par ameter. 



Its value for each parametrization is obtained by fitting it to the H(z) data, 26M30] . 



A. Parametrization I 



As a first parametrization we propose 



.W = -i4M^|. (ID 



2 V<?i + (1 + z)^ 
where to avoid divergences q\ and q2 must be positive-definite. 
Introducing (jlip in ffTUl) . it follows 



3 



"H-M^wT- (12) 

(We note in passing that for q2 = 3 the ACDM behavior is reproduced). By using the 
method outlined in last section in conjunction with the observational data (SN la (557), 
CMB/BAO (7) and H(z) (15)), we fit the three free parameters occurring in (TT2"j) . The 
result is q x = 2.87t{j;£j, q 2 = 3.27 ± 0.55, and H = 70.5±H km/s/Mpc. Table M shows the 
X 2 values of the best fit. 



Data sets 


2 2 
XSN XbAO/C'MB 


Xh 


Xtot 


Xlot/ do f 


Union2+B AO /CMB+Hubble 


542.6 2.6 


17.9 


563.3 


0.96 



Table II. Best fit \ 2 values of parametrization I, Eq. (jlip . The free parameters are qi, q 2 and Hq. 



Figure [T] shows the evolution of the q for the best fit values of parametrization I (solid 
line), with its la confidence region (shadowed area), and the spatially flat ACDM model 

9 



(dashed line) as determined by the WMAP 7-years team 24| (the latter graph is included 
for the sake of comparison), in the interval — 1 < z < 5 (left panel), and the evolution of the 
Hubble function in the interval < z < 3 (right panel). 
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Figure 1. Left panel: deceleration parameter vs. redshift. The shaded area shows the la confidence 
region. Right panel: Hubble function vs. redshift (the observational data are borrowed from Refs. 



2614291] ) . In both panels the solid (red) and the dashed (green) lines are used for the best fit of 



parametrization I (Eq. (fTTl) ) and for the ACDM model with Qmo = 0.27 and Hq = 72.1 km/s/Mpc 



-sec 



24j]-. respectively. The latter graph is shown for comparison. 



Figure |2] depicts the la and 2a contour plots of the pairs (qi, q 2 ) (left panel) and (H , 
go) (right panel). Upon assuming that the expansion is dominated by pressureless matter 
and some other (unspecified) component, non- interacting between each other except gravi- 
tationally, the effective equation of state (EoS) parameter is given by 

2q(z) - 1 

-M- 3(i -»„(,)) • (13) 

Parametrizations of q(z) and w(z) are somewhat equivalent but not quite because in the 
latter case some assumptions about the energy budget of the Universe, as well on the ex- 
istence or not of possible interactions between the different components, have to be made 
while in the former (as in our case) not necessarily. 

Having said this, it is interesting to confront (fT3|) with the widely used Chevallier- 



Polarsky-Linder (CPL) parametrization 3ll. |32| 



w = w + w 1 ——, (14) 
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Figure 2. Left panel: la and 2cr confidence regions of the pair of free parameters (^1,92) of 
parametrization I, Eq. (jlip . Right panel: 1<t and 2a regions of the pair of free parameters (Hq, qo). 
The dot signals the best fit values. 



not far from z = 0, in the redshift range —0.3 < z < 0.3. We restrict ourselves to comparative 
small redshift around z = because it diverges at z — > —1. After numerically linearizing 
our expression for w(z) we get wq = —0.92 ± 0.10 and w\ = 0.3ll ;25> va l ues m ver y good 



agreement with those reported in 
and W\ 



24|, for the CPL parameters, namely: Wq = —0.93 ± 0.12 



-o.38±g : g. 



As Fig. |3] reveals, comparison in the extended interval —0.5 < z < 5 shows that for 
z > 2.5 the evolution of the effective of both EoS disagree in excess of la. This is consistent 
with claims that the CPL parametrization is not appropriate to fit data simultaneously at 
low and high redshifts 



33 



34]. 



B. Parametrization II 

As a second parametrization we propose, 

q ( z ) = _ 3 gi + 1 _ 3( gi + l)« — . (15) 

In this case the Hubble function must be obtained by numerically integrating Eq. (|10p . 
Proceeding as before we obtain, q x = 0.078^;^, q 2 = 0.95lg^, and H = 70.4 ± 1.6 
km/s/Mpc for the three parameters entering H(z). The x 2 best fit values are shown in 
table HU 



11 



0.2 




Figure 3. EoS parameters vs. redshift. The shaded area shows the la confidence region. The solid 
(red) and the dashed (green) lines are used for the best fit of parametrization I (Eq. (jlip ) and for 
the CPL parameterization, Eq. (I14p . both with £Imo = 0.27±0.03. For the CPL parameterization 



the values wq = —0.93 ± 0.12 and W\ 



-0.38" 



-0.66 
-0.65 



obtained by Komatsu et al. 



24| were used. 



Left panel of Fig. H] shows the evolution of q for the best fit value of parametrization 
II (with its la confidence region) and the ACDM model obtained by Komatsu et al. 24]. 
The right panel depicts the evolution of the Hubble function versus redshift in the interval 
< z < 3. 



Data sets 


Xsn Xbao/CMB Xtot Xtot/dof 


Union2+CMB /B AO+Hubble 


542.7 2.7 17.7 563.1 0.96 



Table III. Same as Table II but for parametrization II, Eq. (I15p 



Figure shows the la and 2a contour plots of the pair (q%, q 2 ) (left panel), and (H , q ) 
(right panel). Note that go results more degenerate than in the previous parametrization (as 
well as in the next one). This arises because -as direct inspection shows- in the other two 
parametrizations q^ depends on just one free parameter, q ly while in this parametrization it 
depends on both, qi and g 2 . 
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Figure 4. Same as Fig. [I] but for parametrization II, Eq. (fT5 
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Figure 5. Same as Fig. [2] but for parametrization II, Eq. (fl~5j) . 



Considering the effective EoS parameter (fl~3l) . as in the previous section, we obtain after 
linearization wq = —0.97^^ and wi = — 0.15^'™. As Fig. [6] shows, the evolution of the 
said effective EoS and the CPL in an extended redshift interval is similar to the previous 
one but with the la uncertainty interval significantly wider. 



C. Parametrization III 

The previous parametrization presents the inconvenience of a significant uncertainty in 
go since it depends on the two free parameters. The following parametrization 

12-7== -92- 



1 3 qi e q2 ^ - e 



T+5 



!(«) = -;+7- ; — . (16) 
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Figure 6. EoS parameters vs. redshift. The shaded area shows the la confidence region. The solid 
(red) and the dashed (green) lines are used for the best fit of parametrization I (Eq. (jlip ) and 
the CPL parametrization, Eq. (I14p . respectively, both assuming Qmo = 0.27 ± 0.03. For the CPL 



jarametrization the values wq 
24j were used. 



-0.93 ± 0.12 and wi 



-0.38 



-0.66 
-0.65 



obtained by Komatsu et al. 



avoids this as go depends on q 1 only. 

Again, H(z) must be obtained numerically. Then, proceeding as in the two previous 
instances, we obtain q 1 = 0.36^^1, q 2 = 1.57toH, and H = 70.5~t{ 6 km/s/Mpc. The 
X 2 values of the best the fit are indicated in table [TV] Figure [7] shows the evolution of the 



Data sets 


Xsn Xbao/CMB Xtot Xtot/dof 


Union2+CMB /B AO+Hubble 


542.6 1.7 17.9 563.2 0.96 



Table IV. Same as Table II but for parametrization III, Eq. (I16p . 



deceleration parameter for the best fit parametrization (with its la confidence region) and 
the ACDM model as determined the WMAP team 24j (left panel), and the evolution of the 
Hubble function (right panel). 

Fig. [8] shows the la and 2a contour plots of the pairs (qi, q 2 ) (left panel) and (H , q ) 
(right panel). 
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Figure 7. Same as Fig. [I] but for parametrization III, Eq. (|16p . 
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Figure 8. Same as Fig. [2] but for parametrization III, Eq. (|16p . 

As for the effective EoS parameter ffTB"]) . proceeding as in the previous subsections, we 
obtain wq = —1.01 ± 0.06 and W\ = 0.03 ± 0.16. Again, the evolution of the said effective 
EoS and the CPL in an extended redshift interval is rather similar to the one in Fig. |3j 
thereby, we do not feel it necessary to show it here. 

D. Discussion 

Figure [9] compares the parametrizations. All three yield rather similar results being really 
close between one another from the statistical standpoint {x 2 /dof = 0.96 for all of them). 
However, parametrization II looks somewhat less favored than the other two because of the 
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noticeably wider la region of q vs. z, as seen in the left panel of Fig. HI 
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Figure 9. Deceleration parameters vs. redshift. Solid (red), long dashed (blue) and short dashed 
(green) lines are for parametrizations I, II and III respectively. The graphs of parameterizations I 
and III practically overlap each other. 



All of them are consistent with a present stage of accelerated expansion, never to end 
or slow down. Further, the best fit plots of H(z) and q(z) are alike the corresponding 



plots of the ACDM model as determined by Komatsu et al. 24j . From Table |V] we learn 
that all Hq best fit values are within la of each other and consistent with the Ho value 
reported in [24j. The same holds true for the best fits of the age of the Universe (given 
by to = J °° dz/[(l + z) H(z)]), go and the redshift, z t , at which the transition deceleration- 
acceleration occurred (i.e., q(z t ) = 0), though the central values of the latter are not so close 
between each other as the corresponding values of the other two parameters. At any rate, 
the three of them are consistent with the z t w 0.5 value obtained by Wu et al. using the 



history of the strong energy condition [35j, as well as with the findings of Riess et al. 36], 



Cunha and Lima |5j, and Lu et al. 37] 



As tablelVlshows. the values predicted for Hubble's constant, Hq, by the three parametriza- 
tions are within la between one another and with the value predicted by the ACDM model 
that best fit identical sets of observational data. This is also true for t , q and z t . 
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Param. I 


Param. II 


Param. Ill 


ACDM 


#0 


70-5±i:§ 


70.4 ± 1.6 




70.2 ± 1.4 


to 


13.6 ±0.5 


13.7 ±0.4 


13.6 ±0.2 


13.4 ±0.1 


Qo 


-n 61+ 006 

u - Di -0.07 


U.OD_ 22 


-0.60 ± 0.06 


-0.60 ± 0.03 


Zt 


71+ 014 
u -' 1 -0.17 


n 77 +0.52 
U -' '-0.57 


72+ - 27 

u -' z -0.21 


0.76 ± 0.05 



Table V. Hubble's constant, Hq (in km/s/Mpc), the age of the Universe, to (in Gyr), the deceler- 
ation parameter, qo, and the redshift, zt, of the transition deceleration-acceleration for the three 
parametrizations, and the flat ACDM model fitted to the same data sets. 

A direct and model-independent determination of q(z) in the redshift interval < z < 1 



was carried out by Daly and coworkers 

q(z) = -1 - (1 + z 



who applied the expression [39 1 
d 2 y/dz 2 tt ko ydy/dz 



(17) 



dy/dz 1 + tt k0 y 2 

io the 192 SN I data points of Davis et al. |4fj| and 30 radiogalaxy data points of Daly et al. 
4l| -see Fig. 10 in [38j]. In (JT7J) y(z) = Ho(aor) is the dimensionless coordinate distance, 
r the radial coordinate of the FLRW metric, and Q k o — —k/{H a ) 2 . Notice that it just 
assumes the FLRW metric; i.e., it holds regardless the energy components of the Universe 
or the specific theory of gravity adopted. 

Figure [10] shows the best fit graphs of q(z) of each parametrization superimposed to 
experimental results of Daly et al. j38|. As it is seen, these lines fall within the la region of 
q(z) as determined in Ref. I 



One may ask as to what extent the imposing of the thermodynamic-based, far future, 
constraint q(z = —1) = —1 biases the parametrizations toward the ACDM model. We 
have studied this by letting the value of q(z = —1) as an additional free-parameter and 
fitting it using solely the observational data. The results are: q(z = —1) 
q(z = -1) = 



-0.65 



1-0.21 
-0.50' 



l.lt°7, and q(z 



-0.82i '4 for parametrizations I, II and III, 



respectively. Except for the second one, the quintessence cold dark matter (CDM) model is 
somewhat preferred but, in all the cases, the physically motivated choice q(z = —1) = —1 
results compatible within la. The drawback of letting q(z = —1) free, aside from violating 
thermodynamics, is that the other two free-parameters present a wide degeneracy. 

Likewise, the derived values for the Hubble constant (first row in Table |V]) differ from the 
recently obtained by Riess et al., H = 74.2 ± 3.6 km/s/Mpc, [26], who used 240 Cepheids 
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Figure 10. Top panel: observational constraints, at la confidence level, on the deceleration pa- 
rameter versus redshift, obtained by Daly et al. using the combined sample 192 supernovae and 



30 radio galaxies (top panel of Fig. 10 in 3g]); the thick solid line corresponds to the central 
experimental value of q(z). Solid, dashed, and dot-dashed lines correspond to parametrizations I 
(Eq. (fTT]) ). II (Eq. (fT5|)) and III (Eq. (fTB]) ). respectively. Bottom panel: The same as the top panel 
except that the observational constraints on q(z) were obtained solely from the sample of 30 radio 
galaxies (bottom panel of Fig. 10 in [38f] ) . In both panels parametrizations I and III practically 
overlap each other. 



variables at z < 0.1, by about 6% (but they all agree with the latter at la). We have 
considered this by repeating the analysis of subsection IIIC but this time leaving aside the 
mentioned value of Riess et al. The results now are: H = 70.0t\i, H = 70.0lij, and 
H = 70-llj/g Km/s/Mpc for parametrizations I, II and III, respectively. Thus, both sets of 
results are essentially coincident (they differ by less than 1%). They also agree very well with 
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he Hubble constant value observationally derived by Komatsu et al., 70.4 ± 2.5 km/s/Mpc 



24J, using WMAP 7-year data. So, while there is a significant difference between the Hubble 



constant value of Riess et al. and ours, it is not a substantial one; after all, they agree at 
la confidence level. At any rate, the root of the discrepancy may be rightly traced at the 
difference in methods employed. While Riess et al. essentially used astrophysical data, we 
(as well as Komatsu et al.) resorted to cosmological data instead. 



V. CONCLUDING REMARKS 

In this paper we proposed three different two-parameter parametrizations of q(z) valid 
from the matter era (z 3> 1) up to the infinite future (z = —1), modulo H(z) > 0. These 
rest in the following hypotheses: (i) at cosmological scales the Universe is homogeneous 
and isotropic, thereby well described by the FLRW metric; (ii) in the matter dominated 
era q = 0.5 ; (in) at least at late times the entropy of the Universe is dominated by the 
entropy of the apparent horizon. The second an third hypotheses furnish two fixed points 
(at z 3> 1 and z = — 1, respectively), thereby drastically reducing the ample latitude one 
faces in parameterizing q(z). By smoothly interpolating between these two points one can 
obtain useful parametrizations, but with shrunk arbitrariness. 

Except for the existence of a matter dominated era at early times, the parametrizations 
are independent on any specific cosmological model; and, on the other hand, they are flexible 
enough to accommodate many homogeneous and isotropic models. We constrained the free 
parameters with the latest observational data (SN la, BAO, CMB, and H(z)). Accordingly 
if to accommodate a given cosmological model the free parameters, qi and q2, in the three 
parametrizations should take values widely apart from their respective best fits (which are 
consistent within la with the flat ACDM model), we may confidently discard the said model. 

Thermodynamics in spatially flat (k = 0) FLRW universes demands that q(z = —1) = 
— 1. This provides us with an additional and very useful fixed point to parametrize q(z) 
in a model independent manner. Note that in the absence of a physically-based guidance 
to figure out the value of the deceleration parameter at z = — 1 one is led to choose some 
or other random value. By contrast, in our case we have taken q(z = —1) = —1 on solid 
thermodynamic grounds. 

Albeit we have considered just the particular set of spatially flat FLRW universes, it is 
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not a big restriction at all. Indeed, recalling that in the case of non-flat metrics the area of 
the apparent horizon is given by Eq. (TjQ) it follows that A' = A 2 / (27ca) [H 2 (l + q) + k a~ 2 ]. 
When a — )■ oo and k = —1, the last term on the right hand side is necessarily subdominant 
otherwise one would have A' < 0, contrary to the second law. Hence the condition A' > 
in that limit reduces to the one in flat space, namely, q > — 1. In the positively curved 
case, and again in the same limit, one can assume that H oc a n , n being some real number. 
As it can be straightforwardly checked, the aforementioned last term results, once more, 
subdominant provided that n > — 1 which is the case of most realistic cosmologies. Thus, 
(z = —1, q = —1) is an asymptotic fixed point not only for spatially flat universes but also 
for open universes and for a rather ample set of closed universes. 

Our results suggest that from the era of matter domination onward q decreases monotonously 
with expansion (i.e., dq/dz > 0), the transition deceleration acceleration occurred at a red- 
shift of about 0.7, and that go — —0.6. They do not support recent claims that the cosmic 



34L m-m- 



46] , as well 



expansion is today reverting to a decelerated phase (i.e., that dq/dz\n < 0) 
On the contrary, they show overall consistency with the findings of {3]-[iq|], 3| 
as with those of Serra et al., 47(|. The latter authors showed that the equation of state of 
dark energy has not varied noticeably in the redshift interval < z < 1. 
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